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ABSTRACT 

We follow the time dependent thermal evolution of a white dwarf (WD) undergoing sudden accretion in 
a dwarf nova outburst, using both simulations and analytic estimates. The post-outburst lightcurve clearly 
separates into early times when the WD flux is high, and late times when the flux is near the quiescent level. 
The break between these two regimes, occurring at a time of order the outburst duration, corresponds to a 
thermal diffusion wave reaching the base of the freshly accreted layer. Our principal result is that long after 
the outburst, the fractional flux perturbation about the quiescent flux decays as a power law with time (and not 
as an exponential). We use this result to construct a simple htting formula that yields estimates for both the 
quiescent flux and the accreted column, i.e. the total accreted mass divided by WD surface area. The WD 
mass is not well constrained by the late time lightcurve alone, but it can be inferred if the accreted mass is 
known from observations. We compare our work with the well-studied outburst of WZ Sge, finding that the 
cooling is well described by our model, giving an effective temperature Teff = 14,500 K and accreted column 
Ay ss 10® g cm“^, in agreement with the modeling of Godon et al. To reconcile this accreted column with the 
accreted mass inferred from the bolometric accretion luminosity, a large WD mass > 1. IMq is needed. Our 
power law result is a valuable tool for making quick estimates of the outburst properties. We show that fitting 
the late time lightcurve with this formula yields a predicted column within 20 % of that estimated from our full 
numerical calculations. 

Subject headings: accretion, accretion disks — novae, cataclysmic variables — stars; individual (WZ Sagittae) 
— white dwarfs 


1. INTRODUCTION 


Dwarf novae (DNe) are cataclysmic variables (CVs; 
Warner 1995) that undergo dramatic accretion events in which 
a large fraction of the accretion disk is dumped onto the 
white dwarf (WD) surface. These last for ~ 2 - 20 days 
at accretion rates of ~ 1 O“*^M 0 yr“' and are separated by 
quiescent intervals of ~ 10 days to tens of years. Fol¬ 
lowing D N outbursts , CVs show a decreasing ultraviolet 


flux (e.R.. iLons et alJ 11994 

Gansicke & BeuermannI 119961 

ISion et alJll996HSzkody et al. 

ll998llChengetal.ll2000l). The 


accretion rate then is low enough that the WD surface can 
be directly seen, and spectroscopic observations indicate 
this flux originates from a WD cooling in response to the 
outburst_(|Mate(j&^zko^|1^8^|K ir)linger. Sion, (S^zkodvl 
fT^ iLong et all Il993t ISionI 11993^ ! The flux converges 
asymptotically to a quiescent level set by th e time-averaged 
accretion rate llTownslev & Bildstenl l2005l) . Fq = trsBTrffo^ 


where ctsb is the Stefan-Boltzmann constant and Teff.o is the 
quiescent WD effective temperature. 

The lightcurve of the cooling WD depends on the 
mechanism for heating, which has previously been at¬ 
tribut ed to boundary lay er irrad iati on ilRegevll1985l FPringle 


iRegev 


1995 


1995 


/ & Shara Il989t iGodon. Regev. & Shayiv 
iPophamI 1 1997), compressional heating llSior 
iGodon & SionI l2002h . and shear mixing luminos- 


ity_j| 2 E 2 ^S.£L^Ji 222 li' ^ coyering fraction < 0.1 

llPiro & Bildstenl 1200411 we show in the Appendix that 
boundary layer heating (eyen yia shear mixing) and surface 
adyection haye a negligible impact on the late time cooling 
lightcurye as long as the imposed surface temperatures 
are <10® K. Temperatures higher than this are ruled 
out according to theoretical work on the bou ndary layer 
llPopham & Naravanll 99.51 IPiro & Bildstenll2004ll . We there¬ 
fore focus solely on compressional heating to understand the 
cooling flux released on long timescales. We show that this 
late time cooling is a powerful tool for studying DN outbursts 
as it is independent of the time dependent accretion rate dur¬ 
ing the outburst. Furthermore, WD photospheric temperature 
estimates are less uncertain well after the outbur st since any 
residu al screening material has dissipated (see iLong et 

l2(inl . 

We present numerical calculations following the compres¬ 
sional heating and subsequent cooling of the WD surface in 
§2. We inyestigate the late time cooling analytically in §3, 
showing that it obeys a power law, and not an exponential 
decay (as has been assumed i n other studies, for example in 
iGansicke & BeuermannlfT99^ . Furthermore, this power law 
is not in temperature, but rather in fractional temperature 
perturbation about the quiescent WD temperature, ^Teff/Teff^o. 
We construct a fitting formula (eq. 1201 ') that can be applied to 
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DN cooling observations to constrain both the quiescent flux, 
Fq = csfiTrff Q, and the column of material accreted during the 
DN outburst. Ay = AM/(47rR^) (using eq. 1221 ). where AM 
is the total amount of mass accreted. This late time lightcurve 
depends weakly on the surface WD gravity, and therefore 
cannot directly constrain the WD mass. However, if AM is 
known from other observations, the WD radius (and mass) 
can be inferred. We compare our model to the 2001 July out¬ 
burst of WZ Sge in §4, and show that its cooling is consis¬ 
tent with our model. Our analytic formula fits the late time 
lightcurve and yields an accreted column to within 20% of 
that found from our full numerical calculation. We conclude 
in §5 by summarizing our results and discussing the impor¬ 
tance of further, multi-epoch observations of WD cooling fol¬ 
lowing DN outbursts. An appendix investigates the impact of 
surface heating on the late time cooling lightcurve. 

2. TIME EVOLUTION DURING THE DWARE NOVA OUTBURST 

In this section we present time dependent numerical mod¬ 
els of WD thermal evolution during and after a DN outburst. 
These highlight the importance of deep compressional heat¬ 
ing, which will be useful for constructing the analytic late 
time cooling curves. 


2.1. Quiescent White Dwarf Model 


We begin by describing the quiescent background model 
of the WD envelope, as expected between outbursts when 
the accretion rate is low (M < IO'^Mq yr“'). We focus on 
shallow, radiative surface layers of the WD, which we model 
with plane parallel geometry and constant gravitational ac¬ 
celeration, g = GM/R^. It is convenient to use column depth 
y = -f pdr as a vertical coordinate, where r is the spherical ra¬ 
dius. Hydrostatic balance then becomes F = gy. Throughout 
this paper we give numerical estimates using solar compo¬ 
sition, ideal, nondegenerate gas and Kramer’s opacity. We 
made comparis ons with a more detaile d calculation using 
OPAL opacities (llglesias & Rogersil996h and did not And sig¬ 
nificant deviations from our results. For solar composition 
and gravity g = 10^ cm s“^, surface convection zones are con¬ 
fined to columns y < 1 g cm“^ for Tgff > 11,000 K, and can be 
ignored. 

In quiescence, the thermal profile is set by the underlying 
flux from the core, determined by the time-averaged accretion 
rate (iTownslev & Bildstenll2003li . The WD surface layers are 
described by the radiative flux equation 

3k dy’ ^ ^ 

where F is the outward directed flux, ass is the Stefan- 
Boltzmann constant and k is the opacity, which we assume 
satisfies a power law k = Kop‘‘T^. For Kramer’s opacity a= 1, 
b = -l /2, and kq = 6.5 x 10^^ (in cgs units). Denoting the qui¬ 
escent flux by Fq = assT^f^Q, we integrate equation o to And 
the temperature profile in quiescence 


To(y) = 


(n-l-1)— 

a 


1 l/(4+a-i) 




!/(«+!) 


= 5.7xlo5KgfVo"/^"yf\ (2) 


where gg = g/10* cm s F ()\2 = Tq/IO'^ ergs * cm 
ys = y/10^ gcm“^, n = {3-b)/{\ + a) is the poly¬ 
trope index (n = 13/4 for Kramer’s opacity), and a = 
(l6asB/3KQ){kb/pnipgf = 6.2 x lO'^'^gg* (in cgs units). 


A property of the envelope that is important for understand¬ 
ing the cooling is the local thermal time at a column depth y. 


fth.oCv)^ 


yCpTpiy) 

Fq 


(n+1) 


Fq 




1 /{4+a-b) 


y 


{n+2)/{n+l) 


r, no 2/17^^-15/17 21/17 
= 0.73yrg8' ^o,i2 ^5 > 


(3) 


where Cp ~ Sk^ / (2pmp) is the specific heat, k^ is Boltzmann’s 
constant, and p. = 0.62 is the mean molecular weight in the 
nondegenerate plasma of solar composition. 


2.2. Compressional Fleating During the Outburst 

The entropy equation describing time evolution of the WD 
temperature pro file during the DN outburst can be written as 
dBildstenl 1998ft 


dF 

dy 


~ (2 n 


dT riiT „ ' 

-^ +-(V-Vad) 

ot y 


(4) 


where m = M/(AttR^) is the accretion rate per unit area, 
V = d\t\T/d\t\y is found from equation Q, and Vad = 
(OlnT/c)lnP)ad « 2/5 sets the adiabatic profile. The sec¬ 
ond term on the right hand side arises from the advection of 
entropy Vrdsjdr = {-m/p)ds/dr = mdsjdy by the accretion 
flow, where Vr is the advection velocity and s is the specific 
entropy. 

We solve equations o and 0 implicitly in time using 
backward time differencing for stability and typically 128 to 
512 grid points. The time step is chosen so that the average 
temperature at each grid point changes by a fractional amount 
< 10“^ from one time step to the next. We repeated a number 
of runs with a fractional change of 10“"^ and the results did 
not change. The boundary condition at the base of the layer is 
taken to be constant flux, Fq, and set at a depth (« 10^ g cm“^) 
where the local thermal time is longer than any timescale of 
interest (i.e. longer than the time we wish to follow the cool¬ 
ing). The flux into the core does change due to each DN out¬ 
burst, but by a tiny amount of order the ratio of the mass ac¬ 
creted in one DN outburst to the envelope mass. As we do not 
model the photosphere, the surface boundary condition is set 
at a shallow depth with a sufficiently short thermal time for 
the temperature profile to be constant flux, T oc y'/(«+i) xhe 
initial condition for the outburst is the constant flux profile 
in equation 0. This is evolved forward for a time fdno, with 
mass accretion set by a given function m{t). 

In Figure n we plot the heating profiles during a DN out¬ 
burst at 10 different times, spaced logarithmically from 0.52 
to 52 days. The WD has Teff o = 14,500 K and g = 10^ cm s“^ 
(as appropriate for an M = Q.6Mq WD), and accretes at a 
constant rate m = 0.30 g cm“^ s“' [a global accretion rate of 
6 X 10“*Mq yr“'(R/10^ cm)^], giving a total accreted column 
Ay = 1.3 X 10® g cm“^. We focus on the fractional tempera¬ 
ture and flux profiles, 5T(yf)/TQ(y) = (T(y,t)-To(y))/To(y) 
and SF(y, t)/F q = {F(y, t) - Fq)/Fq, respectively, as we later 
show that these quantities obey power laws at late times. 

The bottom panel in Figure ^ shows the flux perturbation 
as a function of depth, increasing in time from bottom to 
top. Two distinct regimes are seen for 6F/Fq, separated by 
the base of the accreted layer, shown by dots. Above this 
depth, the profiles are determined by the steady state solutions 
to equation 0, while below the perturbations are small and 
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Fig. 1.— Temperature, T, fractional temperature perturbation, ST/To, 
and fractional flux perturbation, |5F/fo|, during accretion with m = 
0.30 g cm“^ s“*, fjno = 52 days, and Teff o = 14,500 K. The 10 curves are 
logarithmically spaced in time from 0.52 to 52 days, from bottom to top. The 
dots in the bottom panel mark the base of the accreted layer, mt. At some 
depths F(y,t) < Fq, so that SF/Fq < 0. Since we plot the quantities loga¬ 
rithmically to emphasize the power laws, we must take the absolute value of 
SF/Fq, which creates the jumps seen in the bottom panel whenever SF/Fq 
passes through zero. 


largely due to adiabatic compression. The middle panel shows 
the fractional temperature perturbation. Near the surface, this 
quantity is fairly constant with a slight decrease with depth, 
but below the base of the accreted layer it falls off as a power 
law ST /To oc y“' for adiabatic compression. Finally, the ab¬ 
solute temperature is shown in the top panel. The flattening 
of the temperature profile near the base of the accreted layer 
is apparent, especially at late times. If the outburst had lasted 
longer, or the accretion rate was higher, a temperature inver¬ 
sion would have developed just below the base of the accreted 
layer. Since t± is short in these shallow layers, such an inver¬ 
sion “smooths out” very quickly once the envelope begins to 
cool and has a negligible effect on the cooling lightcurve. 

The final profile after the DN event can be understood by 
comparing the timescales of equation 0 , which are the ac¬ 
cretion time face = y/'w, the thermal time f* = yCpT/F, and the 
elapsed time f. Above the base of the accreted layer, y <C mt, 
there is a hierarchy of timescales t± ^ face t, so that we take 
d/dt kQ. The flux near the surface is nearly constant, with a 
small decrease with depth due to compressional heating. De¬ 
noting the (constant) surface flux as Fs, we hnd the temper¬ 
ature prohle associated with this, which we denote Ts(y), by 
using equation 0 with Fq replaced by Fg. Substituting Ts(y) 
into equation 0 and integrating we find 

(5) 

This shows how the flux decreases near the surface, as can 
be seen in Figure [2 The true temperature profile is there¬ 
fore smaller than the constant flux prohle, 7j(y), and can be 
estimated by substituting equation 0 into equation 0 and 


integrating, 

T(y)^Tg(y) 



1 

5 + a — b 



mCpTgjy) 

Fs 


(6) 


Setting T^ ~ 0 at the base of the accreted layer, y « ffifdno (see 
Figure 0, we can estimate the Hux Fs as a function of fdno- 
This implies an effective temperature during accretion 



«3.0x lO^Kgg'^^® 


fn \ ^ / fdno 

0.3 g cm“^ s“* / \ 1 day 

which agrees with our numerical models. This power law in- 
dex of 1/15 » 0.067 is similar to the power laws found by 
iGodon & SionI <20031) during heating. We also agree that the 
power law during heating is for the actual effective tempera¬ 
ture (as opposed to the fractional effective temperature as we 
hnd for the cooling). Since this evolution takes place coinci¬ 
dent with a high accretion rate, it will not be observable so we 
do not elaborate on its effects any more here. 

We now derive the prohle for ST /Tq deep in the envelope 
due to adiabatic compression. This is the hrst step in under¬ 
standing late time cooling. Well below the base of the ac¬ 
creted layer, y 3> Ay, the timescales satisfy f ^ face ^ hh- The 
initial temperature prohle from equation 0 is compressed 
adiabatically, giving 




T(y) = To(y-Ay) 



Vad 


^To(y) 


Ay 

l + (Vad-V)^ 

y J 


( 8 ) 


or, a fractional temperature perturbation 


ST 


'(Vad-V) 


Ay 


(9) 


which explains the scaling ST /Tq oc y ^ in Figure 0 This ex¬ 
tends from deep in the envelope up to y Ri Ay = mfdno- 


2.3. Cooling Following the Outburst 


We evolve the hnal prohle from Figure0forward in time to 
see how the WD cools after accretion has halted. Again we 
solve equations 0 and 0 , but now with m = 0. In Figure 
0we show prohles during cooling spaced logarithmically in 
time from 10 to 1000 days after the end of the outburst. Near 
the surface, the prohles show constant perturbations SF/Fq 
and ST /Tq, while at large depths the prohle is unevolved. Just 
as during heating, SF /Fq is an order of magnitude larger than 
ST /Tq, which we explain in ^when we derive self-similar 
solutions for the cooling phase. 

From this example we see that cooling involves a thermal 
wave moving into the envelope, propagating down to a col¬ 
umn y only at a time, f ss fth(y)- This is simplest to show for 
the deep, adiabatically compressed material. Since the pertur¬ 
bations at these depths are small, we approximate f* Ra fth_o. 
Setting fth,o = f we invert equation 0 to hnd the column that 
is just beginning to cool as a function of time. 


yth(0 ^ 


/Pq \ («+l)/(n+2) 

a 

VcJ 

(n+l)Fo 


1 l/(n+2) 


T«+l)/(n+2) 


a.3 X 10^gcm“^ 


-2/21pl5/21 
8 -^ 0,12 



17/21 

, ( 10 ) 






















4 


PIRO, ARRAS, & BILDSTEN 



y[g cm 2] 


Fig. 2.— Profiles for T, ST /To, and \SF/Fo\ during the cooling following 
the outburst presented in Figure^ in which a column Ay = 1.3 X 10*’ g cm ^ 
was accreted in /duo = 52 days. The 10 curves are logarithmically spaced in 
time from 10 to 1000 days, from top to bottom. 


where t is measured from the end of the outburst. Substituting 
this result fory in equation 0 , we find that the late time cool¬ 
ing should evolve according to ST/Tq ~ (Vad-V)Ay/yth(0 oc 
f-(n+i)/(n+2) .^g shall see in ^ this estimate is correct up 

to a numerical factor that can be computed. 

In the Appendix we consider the effects of boundary layer 
heating on this late time cooling. We find that surface heating 
only propagates a finite distance into the envelope, while com- 
pressional heating causes a ST/Tq oc y“' perturbation down to 
arbitrarily large depths, so that at very late times compres- 
sional heating always dominates. In practice, this may occur 
so long after the outburst that during the observational pe¬ 
riod both heating mechanisms must be considered for proper 
modeling. This happens when the boundary layer is too hot, 
so that a diffusion wave propagating from the surface does not 
have sufficient time to reach the background quiescent profile 
during he ating , resulting in a critical boundary layer tempera¬ 
ture (eq. IA4I '). 


TbI.t 


10'’K 


2/17„2/17 
?8 ^ 0,12 


Ay 


10® g cm 2 


4/17 


( 11 ) 


This critical temper ature is close to, but most likely ruled 


out by observations i 

Mauchel2004h and theoretical modeling 

dPonham & Naravan 

1995HPiro & BildsteiJ2004. 


3. SELF-SIMILAR COOLING 


In 32.31 we found that the fractional temperature perturba¬ 
tion decreases as a simple power law with time for late time 
cooling after an outburst. We now provide a rigorous deriva¬ 
tion of the power law index and calculate the prefactor. 

The discussion in §2.3 suggests that the late time cooling is 
governed by the self-similarity variable 




fth,o(y) 

1/2 

y 

t 


.yth(f). 


1 («+2)/2(n+l) 


( 12 ) 


where t is measured from the end of the outburst, and we have 
introduced the 1/2-power for calculational convenience. At 
large depths, the initial temperature profile is given by equa¬ 
tion (|3, so we make the ansatz 


(57’(y,f) = ro(y)(V,d-V)^/(0, 

y 


(13) 


where /(0 is an undetermined function that contains the time 
dependence of the temperature perturbation. To understand 
the late time cooling we now find / for large t (i.e. in the limit 
C-0). 

Taking perturbations of equation 0 about the quiescent 
profile (with m = 0) gives 

dST dSF 
dt dy ’ 

where the perturbed flux is 
d 


SF 

— =(n-l-l) . 

Fq dlny \ Tq 


(ST\ ST 

I — ) +(4+a-fe)—. 


(15) 


We substitute the perturbed temperature from equation O 
into equations d and d, which are then combined to give 
a single equation for /, 






d^[2 4(n+l) 2(n+l) 


(16) 


4(n-l-l)^ 

In the limit ^ ^ cx), / ^ 1 by definition. In the other limit 
of ^ ^ 0, equation ( fT^ has a finite solution / = jy^ 2 (n+i)/(n+ 2 )^ 
where p is a constant of proportionality determined by the full 
numerical solution. Since equation d contains unwanted 
solutions that are strongly divergent, we solve it as a matrix 
equation, finding rj = 0.141 for Kramer’s opacity. The late 
time temperature perturbation is then 


ST(y,t) Ay 

^i7(Vad-V)- 


(17) 


Toiy) yth(0 

which is what we guessed in §2.3, up to the factor p. 

To make comparisons with the observed late time cooling, 
we relate the above analysis to the fractional Teff perturbation, 
STeff/Tefffl = (Teff- 7’eff_o)/7’eff,o « SF/ 4Fq. At late times, 

-l^2(«+l)/(«+2) _ jjj (18) 


Tq 


(xy 


Equation d then relates the temperature and flux perturba¬ 
tions 

SF ST 

^{4 + a-b) — , (19) 

Fq Tq 

where the factor A-+a-b = ^.5 explains the order of magnitude 
difference between SF /Fq and ST /Tq. Both SF /Fq and ST /Tq 
are constant in space in the limit t fth(y)- We put these 
results together to express the change of ST^a in the linear 
limit in the form 

Teff.O^V t ) 


( 20 ) 


where t is the time since the outburst ended, the power law 
indexis (n-l-l)/(n-l-2)= 17/21 «0.81, and the characteristic 
late time cooling timescale is 

- (n+2)/(«+l) 


A-\-a-b _. 

^late ~ ^ ^ (vad~ v) 

i 0.024 fth,o(Ay). 


fth,o(Ay) 


( 21 ) 
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Observations that fit fiate then constrain 
Ay = 1.8 X g cm“^ 

17/21 / rj. 


X8s 


2/21 ( hate 

1 day 


VlO^K 


60/21 


( 22 ) 


the total accreted column during the outburst. 

Equations and Mil are our central results. We make 
the following remarks on their use for understanding observa¬ 
tions: 


• Observations should be plotted as in order 

to test whether a power-law is observed. Equation (|20j 
has three parameters, Teff o, hate, and n, the polytrope 
index. If the wrong Teff o is used, the late time cool¬ 
ing does not look like a power law with index Ri 0.81, 
giving a useful constraint on Teff o. Eigure|3 shows the 
result for (reff(f)-2i:7’eff^o)/2i:Teff,o from the full cooling 
code, where x = 0.9,1.0,1.1. At late times, the analytic 
formula given by equations (|20j and (I22> agrees well. 
When the wrong quiescent effective temperature is used 
the lightcurve either does not look like a power law (as 
in the case of too large of a Teff.o, i.e. x = 1.1), or shows 
a power law that differs considerably from a power law 
index of 0.81 (as in the case of too small of a Teff o, i.e. 
x = 0.9). 

• If the temperature perturbation ^Teff/Teff o ^0.1, effects 
second order in (^Teff/Teffp cause the lightcurve to devi¬ 
ate slightly from a power law, as is apparent in Eigure|3 
for X = 1.0 and t < 400 days. However, adjusting Teff^o 
until the data is nearly a power law with index 0.81 and 
applying the analytic formula gives quick results accu¬ 
rate to ~ 20%. We discuss this point in more detail in 
€1 

• Given a measurement of hate, we can constrain 

using equation (|23. If we allow WDs in the mass range 
0.6-1 .2 Mq, the factor can vary by ~ 6^/^' « 1.2. 
Therefore, if the WD mass is not known, we can con¬ 
strain Ay, the accreted column to 20%! 

If multiple outbursts are observed from the same system, 
the lightcurves can be compared to infer the relative amount 
of mass accreted. We combine equations (|20j and (|23 to find 

AMi _ ST,s,i(t) 

am 2^ sns^ty ^ 

where AM is the amount of mass accreted in an outburst 
and the subscripts refer to the two outbursts. In this way 
disk instability models can be tested, including understand¬ 
ing the differences between normal DN outbursts and super¬ 
outbursts. In Eigure|4]we show how this could be done by 
comparing the lightcurves of two outbursts, one accreting at 
a rate of 0.30 g cm“^ s“' for fdno = 52 days and the other at 
a rate of 0.15 g cm“^ s“' for fdno = 26 days, but both with 
Teff,o = 14,500 K, so that four times the mass is accreted in 
one outburst in comparison with the other. In the top panel we 
compare the two resulting lightcurves, and in the bottom panel 
we look at the ratio of these two lightcurves, which shows how 
this ratio asymptotes to the accreted mass ratio at late times. 
Eor this to work correctly, the two lightcurves must be cor¬ 
rectly positioned in time, with f = 0 corresponding to the end 
of the outburst. Eortunately, at late times (when this works 
best) it should not be difficult to do this with enough accuracy 
(within ~ 10 days) for useful results. 



time[days] 


Fig. 3.— The Tgff perturbation during cooling from the profiles shown 
in Figure 1^ The three solid lines labeled x = 0.9,1.0,1.1 represent (reff(t) — 
■^7'eff^o)/(-'^7’eff,o), to show how a power law with index 0.81 only results if 
the correct Tgff o is used. The analytic formula given in equation I20i is shown 
by the dashed fine. 



Fig. 4.— A comparison of two outbursts, one accreting at a rate of 
0.30 g cm“^ s“^ for /dno = 52 days and the other at a rate of 0.15 g cm“^ s~^ 
for ?diio = 26 days, but both with Tgff o = 14,500 K. This results in four times 
the mass being accreted in the former outburst in comparison with the lat¬ 
ter. In the top panel we show the two resulting ST^ff/Tsfifi lightcurves, while 
in the bottom panel we plot the ratio of these two lightcurves. This shows 
that at late times the ratio asymptotes toward a value equal to the ratio of to¬ 
tal accreted masses. A similar study would be useful when considering the 
multiple outbursts from a given CV. 
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Fig. 5.— Contours of constant (dashed lines) from fitting the eight tem¬ 
perature measurements during the cooling of the 2001 July outburst of WZ 
Sge. We vary both tfi and Teff o over a grid of our numerical models, but keep 
tdno = 52 days fixed. The star marks our best fit model, m = 0.30 g cm“^ s“' 
with Teff.o = 14,500 K. 


4. THE 2001JULY OUTBURST OF WZ SGE 

The most detailed observation of a decay after a DN 
outburst is from th e 2001 July 23 outburst of WZ Sge 
Jlshioka et alJ 120011) . providing an excellent test of our 
cooling model. The outburst lightcurve is complicated 
JPatterson et alJ200^ . initially exhibiting a period of high M 
for 24 days (with a steady decline), a sharp drop in M for 3 
days, and finally quasi-periodic accretion from days 29 to 52. 

The surface temperature of t he WD during cooling has been 
estimated three different ways (iLong et al.l2004tlGodon et al.l 
l2004j) . For the sake of comparison, we fit the median tem per- 
ature measurement (denoted “7i” in lGodon"^ ani2004h . and 
assume an error 1000 K (Gansicke 2004, private commu¬ 
nication). To constrain m and Teff.o we calculate using the 
eight temperature measurements of the cooling over a grid of 
our numerical models, spaced in intervals of 0.05 g cm“^ s“' 
in m and 500 K in T^ff o- All models assume fdno = 52 days 
with the accretion rate constant during this time. In Fig¬ 
ure 13 we show contours of constant x^^ which favors ni = 
0.30 g cm“^ s“' and Teff^o = 14,500 K (shown by a star). 

In Figure |3 we plot the numerical lightcurve of this fa¬ 
vored model along with the corresponding analytic fit using 
hate = 56 days. At late times the measurement errors are of or¬ 
der the temperature perturbations, which at face value should 
make comparisons to our numerical models difficult. Fortu¬ 
nately, the data show a clear trend that closely follows the 
general features of our numerical lightcurve, which strongly 
suggests that the cooling is due to compressional heating. It 
has a shallow slope at early times and then becomes a power 
law at late times with Sr^ff/Teff o cx f“*^-** as we predict. Our fit 

implies an accreted column Ay = 1.3 x 10^ g cor¬ 

responding to an average accretion rate during the outburst 



time[days] 


Fig. 6. — The fractional temperature during the 2001 July outburst of WZ 
Sge in comparison the cooling of a compressionally he ated envelope. The 
circles show temperature ‘T/,” from IGodon et all 1200411 . with the error bars 
of 1000 K. The solid curve is our numerical calculation with Tgff q = 14,500 K 
and constant m = 0.30 g cm“^ s“* for 52 days, giving Ay = 1.3 x 10^ cm“^. 
This con'esponds to an average accretion rate of (M) ^ 10 ~^Mq yr~^ for an 
M = I.IMq, /? = 5 X 10^ cm WD. The analytic formula given in equation 
for hate = 56 days is shown by the dashed line. 

(M) w 10“* Mq yr~^ (assuming a WD radius of 5 x 10^ cm, 
or mass M = 1.1 in reason able agreement with the M’s 
estimated by iFong et ^ J2003h of (1 - 3) x lO"*^ Mq yr“' 
(measured when the outburst was especially bright). We 
therefore favor a massive WD (small radius) to get an M 
closer to these measurements. Suc h a mass estimate is c on- 
sistent with t he value estimated b y iPatterson et alJ i2002l) of 
1.0 ± 0.2Mq. iGodon et all lliol also find reff,o = 14,500 K 
and Ay w 1.3 x 10® g cm“^ (estimated from their Fig. 2) are 
needed to explain the cooling curve. Note that they include 
a prescription for boundary layer heating, while we have ig- 
nored this e f fect. The agreement between our results and 
iGodon et alJ ll2004h indicates boundary layer heating is not 
important in the late time lightcurve (also see Appendix). 

If we fit the last five out of eight temperatures measure¬ 
ments using the analytic formula from equation (|20|) instead 
of our numerical calculations, we find the best fit for hate = 
39 days with Teff.o = 15,000 K. For this fit we assume a power 
law index of 0.81 because the current data does not allow us 
to treat it as a free parameter. This index value is fairly robust 
since both our work here and similar calculations using OPAF 
opacities give similar results. Using equation i22t . this fit im¬ 
plies an accreted column Ay = 1.1 x 10® g within 

20% of the full numerical lightcurve value. This shows the 
usefulness of our fitting formula for making quick estimates 
of the outburst properties. In general, equation (El will over¬ 
predict Teff.o because the lightcurve is merely approaches the 
power law asymptotically (as shown in Figure 13, so that its 
predicted columns are lower limits. 

5. DISCUSSION AND CONCLUSIONS 
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We presented an investigation into the cooling of a WD af¬ 
ter a DN outburst, focusing on late times when the luminosity 
from the WD is dominated by cooling of adiabatically com¬ 
pressed deep layers. We find that the lightcurve is a power 
law, not in Teff(t), but instead in ST^ffit)/T^s^u. This results 
in a fitting function, equation (I20> . that can be used to con¬ 
strain Teff^o and the accreted column. Ay (eq. ES), from 
the observations. The sensitivity to Teff o may be useful if 
outbursts are too closely spaced for the WD to ever cool to 
the quiescent temperature. C ombining this with the work of 
iTownslev & BildstenI j2003h . would allow one to get a better 
handle on the long term accretion history over the CV’s life¬ 
time. Previous studies that do this require careful measure¬ 
ments of Teff o, but as our work shows, compressional heating 
during DN outbursts can have a long lasting impact on the flux 
leaving a WD. In our example models we find a 10% devia¬ 
tion from the quiescent flux 1000 days after the outburst has 
ended! 

Our study highlights the importance of multiple measure¬ 
ments for a given outburst, over many epochs so as to tightly 
constrain the late time cooling. Many DN outbursts have only 
two measurements of their cooling lightcurve, so we cannot 
extend our analysis to these other systems. We made compar¬ 
isons with WZ Sge because it is the best opportunity to test 
whether compressional heating is occurring and to see how 
well Ay and Teff o can be constrained, but even in this case the 
measurement errors are frustratingly large. It is encouraging 
that our work compares favorably to the outburst from this 
object, and we look forward to new, more detailed measure¬ 
ments that will test the effectiveness of our results. Multiple 
measurements would be a powerful tool in conjunction with 
our fitting function to study the properties of DN outbursts, 
especially in comparing multiple outbursts from the same sys¬ 
tem using equation J23V 

The scalings and arguments we use may help to under¬ 
stand the heating of stellar surfaces in other environments. 


Such effects are observable whenever fth(Ay) is much longer 
than the characteristic viscous timescale in the outer accre¬ 
tion disk, and symbiotic binaries may be the most interesting 
case among such systems. Symbiotic binaries exhibit a vari¬ 
ety of outburst types, one of which is commonly referred to 
as a “classical symbiotic outburst” (see Sokoloski 2004 for a 
recent review). In these outbursts the optical brightness in¬ 
creases by one to a few magnitudes over weeks or months and 
then decays over a timescale of months to years. The cause 
of classical symbiotic outbursts remains a mystery, but such 
timescales are characteristic of the DN outbursts we study 
here. The application of our work to these systems would 
test whether these outbursts are also accretion events. Sym- 
biotics have a number of subtleties that must be correctly in¬ 
corporated before they can be modeled, since they accrete at 
M ~ 10“^- 1Q~^Mq yr“*, considerably higher rates than the 
DN systems. Consequently, there is a hotter boundary layer 
and a thicker disk, both of which increase the influence of 
surface heating (see Appendix). Another differ ence is that 
the temperature inversion we quickly mention in 32.2l is more 
pronounced when this much mass is accreted, but most likely 
it will still have a negligible effect on the late time cooling. 
The Teff o is higher for symbiotics because of a larger time- 
averaged accretion rate and possible steady nuclear burning 
dSion & Starrfieldl 1994ft . and this helps to minimize any tem¬ 
perature inversion. 

We thank Boris Gansicke for providing comments on a pre¬ 
vious draft of this paper and Danny Steeghs and Dean Towns- 
ley for thoughtful discussions. We also thank the anonymous 
referee for helpful comments and criticisms. This work was 
supported by the National Science Foundation under grants 
PHY99-07949 and AST02-05956, and by the Joint Institute 
for Nuclear Astrophysics through NSF grant PHY02-16783. 
Phil Arras is an NSF AAPF fellow. 


APPENDIX 

SURFACE HEATING 

To this point we have neglected the effects of surface and shear heating. We now include the effect of a hot layer heating the 
WD surface, during the outburst and the subsequent cooling. Numerical runs with and without compressional heating are done 
for comparison, and we present analytical results for the case of surface heating alone. The result of our analysis is that a surface 
temperature > 10® K is needed for an appreciable deviation from the analytic formula in equation (I20> . We conclude that surface 
heating has a negligible effect on late time cooling because the boundary layer is not sufficiently hot. These results also apply to 
heating due to the shearing of the surface accretion flow a gainst the WD surface, because the shearing takes place at a column 
much less than those of interest for compressional heating (fPiro & Bildstenll2004). Another mec hanism that heats the surface is 
the advection of hot boundary layer fluid in the accretion flow as considered by Ponhaml Jl997h . Such heating can only occur 
down to a column of Ay, and furthermore, this fluid quickly cools in the shallow surface layers near the boundary layer (because 
fth is so short there), long before it can even reach this depth. Advection is therefore washed out by the surface heating at the 
depths of int erest, which are many orders of magnitude deeper. 

We follow F^ingld jl988h and model boundary layer irradiation by enforcing a surface temperature 7bi during the outburst at 
the outermost grid point. However, surfac e heating front irradiation only o ccurs at low latitudes u nderneath the boundary layer. 
This has been shown both observationally JMauchel2004h and theoretically jPrro & Bildstenl2004l) to cover a small fraction of the 
WD surface area, Cbi ^ 0.01 -0.1, while compressional heating occurs over the whole star. Hence, we assume that the lightcurve 
is generated by a fraction ey of the star undergoing surface and compressional heating, and a fraction 1 - Cbi with compressional 
heating alone. In more detail, if Teff i includes only compressional heating and Teff ^ includes compressional and surface heating, 
we make a combined lightcurve as 7’eff= [(l-CbOT^f j-l-ebi7’gff2]*^"^- Other possible heating effects, such as irradiation through the 
radial flux in the disk ( lRegevll983h or direct irradiation from the inner region of the accretion disk, will be of the same order as 
the Cbi w e con sider here, and therefore do not need to be considered separately. 

Figure lATl shows the temperature, temperature perturbation, and flux perturbation underneath the boundary layer at the end 
of the outburst, just before cooling begins. As 7bi is increased, the temperature profile is affected to larger depths. However, 
runs including both compressional and surface heating always show the ST/To ocy”' scaling at large depths, so we conclude: 
surface heating can only propagate in a finite distance whereas compressional heating affects the material at arbitrarily large 
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Fig. A7.— Temperature, fractional temperature perturbation, and fractional flux perturbation at the end of a 52 day outburst in which the star is heated by a 
boundary layer and/or compressional heating. The lines are (m[g cm“^ s“*],rbi[K]) = (0.30,0.0) (solid), (0.30,1.0 x 10^) (dotted), (0.30,2.0 x 10^) (short dash), 
(0.0,1.0 X 10^) (long dash), and (0.0,2.0 x 10^) (dot-dashed), all starting with a quiescent background model with T^ff o = 14,500 K. 


depths. Hence, at late times the cooling must eventually be given by the analytic formulas in equations (|20|l and (|23. In practice, 
for sufficiently high 7bi, this asymptotic scaling may occur so long after the outburst that for the entire observation period both 
heating mech anisms still affect the light curve. We will come back to this point. _ 

Figure lA^ shows cooling after the outburst. An accretion rate m = 0.30 g cm“^ s“* was used for all curves, as in Figure lATl The 
curves with surface heating use 10% of the surface area being heated by the boundary layer, a plausible upper limit. The three 
curves represent compressional heating alone (solid line), and compressional and surface heating, weighted by the surface area, 
for 7bi = 1.0 X 10^ K (dotted line) and 7bi = 2.0 x 10® K (dashed line), respectively. Surface heating is quite sensitive to Tbi, as the 
curve with 7bi = 1.0 x 10® K is identical to that with no boundary layer heating, while the curve with 7bi = 2.0 x 10® K shows 1% 
differences in Teff at 1000 days. 

For the hducial accreted column used in this appendix, which is the column inferred for the July 2001 outburst of WZ Sge 
in surface heating does not affect the lightcurve when the surface area is taken into account. As we have probably used an 
overestimate of the covering fraction Cbi, one would have to increase 7bi to much larger values to com pete with compressional 
heating at late times. Such high boundary layer temperatures are inconsistent with models calculated bv fPiro & Rildsten1(l2004l) 
for M « 2 X 1O“*M0 yr“', which is still more than the average accretion rate we infer in §4 for during the outburst of WZ Sge of 
10 “*^Mq yr“'. 

We now explain the qualitative features seen in Figures IATI and lA8l The hot boundary initiates a thermal wave propagating 
into the star. The column depth to which this nonlinear diffusion wave propagates is estimated from equations o and 0 


ydiff(f) 


T “ ) 


1.1 X 10^ g cm ^ 


Zbi \ 

10® k) 


17/6 



(Al) 


where Tbi e = 7bi/10®K. At depths smaller than ydiff(0, the temperature is nearly 7bi and the flux is approximately constant with 
the value 


rA+a-h 




• 10'4 erg cm ^ s ' 


Tbi 

10® K 


1/2 


1 day 


-2/3 


(A2) 
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Fig. A8.— Teff perturbation as a function of time after end of the outburst of duration 52 days. Lightcurves were generated for ini[g cm“^ s“^],7bi[K]) = 
(0.30,0.0) , (0.30,1.0 X 10^), (0.30,2.0 x 10^). These lightcurves were then weighted by the fraction (10%) of surface area undergoing heating, as described in 
the text. The curves represent compressional heating alone (solid line), compressional and surface heating for Tbi = 1.0 x 10^ K (dotted line) and T hi = 2.0 x 10^ K 
(dashed line), all for a quiescent background model with Teff^o = 14,500 K. The dot long-dashed line is the analytic formula from equations I20i and HU. 


Given a sufficient outburst duration, ^ 3 yr j2^(7bi/10® the thermal wave will penetrate to a depth at which 7bi 

is equal to the background temperature profile. This occurs at a column 


ybi = 


(«+!)■ 


^0 


1 -l/(a+l) 


r«+l 


:1.2 xlO'^gcm-2g-'/V-,f ^ 


(A3) 


A simple estimate for when surface heating will affect the late time cooling is to set ybi equal to the accreted column Ay. At 
this point, the temperature perturbation due to compressional heating becomes small (5T /Tq ^ 1), and additional boundary layer 
heating can delay the late time power law. We hnd the critical boundary layer temperature 


Tblr 


'10‘’K 


2/17 j,2/17 
?8 ^ 0,12 


Ay 


106 g( 


4/17 


(A4) 


consistent with Figure IA7I This derivation only considers layers directly below where surface heating is occurring, so that if 
ebi ^ 1 then Tbi crit must be somewhat larger. Nevertheless, this still provides a useful lower limit for when the boundary l ayer 
may be important. It i s interesting that this 7hi nri f is comparable to, but m ost likely ruled out by, observations auc hel2nn4 and 
theoretical modeling dPonham & Naravanll99-'i : IPiro & Bildstenll2004l) of boundary layer temperatures. If 7bi is larger t han t his 
critical value, the cooling prohle does not reach our analytic result in an observable amount of time, as shown in Figure IXSl so 
that if our htting formula is incorrectly applied it will result in an overprediction of Ay. 

Now consider surface heating at large depths, ignoring compressional heating for simplicity. At depths y ^ ybi, the temperature 
perturbations become small, and linear theory is again applicable. Using the same method as in we insert the ansatz ST /To = 
/(0 into equations o and 0 , where ^ = [fth.oCr)/?]*^^ uses the thermal time of the background state. We can derive a self¬ 
similar equation for the temperature perturbation at large depths during heating, which we do not present here. As we desire the 
temperature perturbations deep in the envelope, we take the 1 limit, hnding the simpler equation 


2(n+l) df 
di'^^ (n + iy-^di^ 


(A5) 
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with a solution at f = fdno of 


5T 


/ocC ‘exp(-A^2^ 


^dni 


-I 1/2 


fth,o(y). 


exp 


-A 


fth.oCv) 


^dn 


(A6) 


where A = (n+ l)/(n + 2)^ « 0.154. He nce the temperature perturbations “tunneling" further than yti are exponentially small, as 
seen in the rapid dropoffs in Figure lATl The small factor A <C 1 implies that the dropoff occurs at a column an order of magnitude 
larger than one would naively expect. 
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